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I. INTRODUCTION 

The Landau (or Lorenz) gauge gluon and ghost 
propagators in pure Yang-Mills theories or in full QCD 
have attracted much interest for many years. 

One reason for this interest is that the behavior 
of these propagators in the infrared (IR) region has 
been related to gluon and quark confinement [H, 0, Q . 
In particular, Zwanziger has argued that the gluon 
propagator should vanish in the IR limit [1,0,01 while 
the ghost dressing function should become singular. 

Another reason for this interest is connected with 
the importance of the momentum dependence of both 
propagators, especially in the (deep) IR-region, for 
the phenomenological analysis of experimental data. 
Many years ago Parisi and Petronzio Q have pointed 
out that non-zero effective gluon mass (or dynami- 
cally generated gluon mass) m g is important to re- 
solve some discrepencies in low-energy tests of QCD, 
as, e.g., ratios of widths of J/^p. Since then a num- 
ber of papers has been dedicated to phenomenologi- 
cal studies of such processes as J/ip — > jX, 7 — > ir° 
transition, non-leptonic B meson decays, etc., where 
a non-zero value of the effective gluon mass m g plays 
a crucial role (for an incomplete list of references see, 
e.g., the recent papers @, Q). 

Let us note that in order to obtain a reliable value of 
m g one needs to know the continuum gluon propaga- 



tor D(p) in the deep infrared region. The definition of 
the mass m g is based on the hypothesis of a pole-like 
behavior, i.e. ~ l/(p 2 + m2 g ), of the gluon propaga- 
tor at small momenta. In this case the effective gluon 
mass defines the infrared mass scale miR. 

Gauge- variant QCD Green functions may serve also 
as input to bound state equations as Bethe-Salpeter 
or Faddeev equations for hadron phenomenology [1, 
[Tol [ill ]. Moreover, at large momenta they should al- 
low a determination of phenomcnogically relevant pa- 
rameters such as A-jy-g- or condensates {A 2 ), (tpt/j) 
by fitting lattice data to continuum expressions ob- 
tained from operator product expansion and perturba- 
tion theory [l2], [l3| in a certain MOM scheme [3, [HI ■ 



The search for intertwined asymptotic gluon and 
ghost propagator solutions of Dyson-Schwinger (DS) 
and functional renormalization group (FRG) equa- 
tions showed the existence of infrared solutions ex- 
hibiting a power-like scaling behavior [^, [l6l \vA 

HE 



IT9L I20L 21 . l22l. However, as has been pointed out 
111 there are also regular so-called de- 

coupling solutions providing an IR-finitc limit of both 
the gluon propagator and the ghost dressing function. 
In [27| it has been argued, that it seems to be a ques- 
tion of IR boundary conditions posed on the ghost 
dressing function, what kind of solution one has to 
select. 
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On the lattice, over the last decade extensive stud- 
ies of the Landau gaugegluon and ghost propagators 
have been carried out [H US H Wl [H, [33, O, S 
HE H3, Hl|- In the meantime lattice computations 
have reached lattice volumes even with a linear ex- 
tension O(10fm) in order to discriminate between the 
above mentioned IR solutions. In this way lattice 
QCD has been found to support the decoupling so- 
lution [H, EH EH, S3, El- This is true as long as 
one relies on finite-box periodic boundary conditions 
and on a gauge condition requiring the Landau gauge 
functional to take extrema as close as possible to the 
global extremum. 

The Gribov copy influence still remains a serious 
problem in the lattice calculations, at least, in the 
deep IR-region. We believe that the correct gauge 
condition should require the Landau gauge functional 
to take extrema as close as possible to the global ex- 
tremum . Indeed, 

(i) a consistent non-perturbative gauge fixing pro- 
cedure proposed by Parrinello-Jona-Lasinio and 
Zwanziger (PJLZ-approach) [48|, El presumes 
that the choice of a unique representative of 
the gauge orbit should be through the global ex- 
tremum of the chosen gauge fixing functional; 

(ii) in the case of lattice compact U(l) gauge the- 
ory in the weak coupling (Coulomb) phase some 
of the gauge copies produce a photon propaga- 
tor with a decay behavior inconsistent with the 
expected zero mass behavior [5(| [EH [E2 ■ The 
choice of the global extremum permits to obtain 
the physical - massless - photon propagator. 

For all practical purposes the system of DS and/or 
FRG equations has to be truncated. The details 
of truncation influence the behavior of the Green 
functions especially in the non-perturbative momen- 
tum range around IGeV, where the Landau gauge 
gluon dressing function exhibits a pronounced max- 
imum. Therefore, reliable results from first princi- 
ples to compare with are highly welcome. On the 
lattice, for finite volumes such results can be ob- 
tained and directly compared with finite volume DS 
and FRG results [53j,l54j. To our knowledge, on the 
lattice a systematic continuum limit determination 
of the Landau gluon and ghost propagators for var- 
ious fixed physical volumes in the momentum range 



0.3 GeV 



< 



P 



< 



10 GeV is still missing. Such an 



evaluation has to make sure that Gribov copy effects, 
lattice artifacts and multiplicative renormalizability 
are sufficiently under control. 



Here we present such a study for SU (2) pure gauge 
theory as a continuation of our investigation in [42j . 
For gauge fixing we rely on the Landau gauge with a 
gauge condition requiring the gauge fixing functional 
to take extrema as close as possible to the global ex- 
tremum. For this aim we employ the simulated an- 
nealing algorithm [42l . I55L l56| in combination with 
non-periodic Z{2) gauge transformations ("Z(2)~ 
flips") [13] • With respect to the latter we search 
within all 2 4 — 1 6 global Z(2) Polyakov loop sec- 
tors. We concentrate on the gluon propagator and 
will present data for three different physical volumes 
up to ~ (7.3fm) 4 . As done in a preliminary man- 
ner in [35l | we check for multiplicative renormalizabil- 
ity and provide results for the renormalized propa- 
gator and dressing function, respectively, which can 
be considered already to be continuum ones. In Sec- 
tion m we introduce the observables to be computed. 
In Section [m] some details of the simulation are given, 
whereas in Section lTVl we discuss the effect of improved 
gauge fixing. In Section [V] we present our numerical 
results. Section [VI] is dedicated to the discussion of 
the deep infrared region and to the definition of an al- 
ternative infrared (mass) scale toir. Conclusions will 
be drawn in Section [VIll 



II. THE GLUON PROPAGATOR: 
DEFINITIONS 

In order to generate Monte Carlo ensembles of non- 
gauge-fixed SU(2) gauge field configurations we use 
the standard plaquette Wilson action 



S = p 



EE 

P = 4/5o 



1 - - Tr 
2 



(1) 



go denotes the bare coupling constant, U xfl S SU{2) 
are the link variables. The latter transform under lo- 
cal gauge transformations g x as follows 



U 



XfJ, 



U° = 9tU x ^g x+fi , g x g SU{2) . (2) 



The standard definition [58|] for the dimensionless lat- 
tice gauge vector potential A x+ ^/2 : ^ is 



A 



-(u x 

2i V 



^ X fl 



Ac+/t/2,M 9 ' (3) 



This definition, which is not unique, can influence the 
propagator results in the IR region, where the contin- 
uum limit is hard to control. 

In lattice gauge theory the usual choice of the Lan- 
dau gauge condition is [58| 



1 For recent alternative attempts see Q, [45|, [4(| [47| . 



(dA) x 



4 

E 



(Ax+jX/2-ij, - A x -jX/2;ij) — , (4) 



3 



which is equivalent to finding a local extremum of the 
gauge functional 



Fu(g) 



_L V - Tr U" , 

4V ^ 2 Xfi 



(5) 



with respect to gauge transformations g x . V = L 4 
denotes the lattice volume. The manifold consisting 
of Gribov copies providing local maxima of the func- 
tional ([5|) and a semi-positive Faddeev-Popov opera- 
tor is called the Gribov region while that of the 
global maxima is called the fundamental modular re- 
gion (FMR) A C Vt. Our gauge fixing procedure is 
aimed to approach A. 

The (unrenormalized) gluon propagator D and its 
dressing function Z are then defined (for p ^ 0) by 



-(A;(k)AK-k)) 

So 



Suu 



Pn Vv \ xa b 



P 



5 ab D(p) , 



Z(p) = D(p) p 2 , 



(6) 
(7) 



where A(k) represents the Fourier transform of the 
gauge potentials defined by Eq. (J3]) after having fixed 
the gauge, a denotes the lattice spacing. The physical 
momenta p are given by p^ — (2/a) sin iirk^/L), k^ £ 
(— L/2,L/2]. For p ^ 0, one determines Dip) accord- 
ing to Eq. (0 



a— 1 fi—1 



(8) 



whereas the "zero momentum propagator" Dip = 0) 
is defined as 



3 4 



(9) 



a— 1 /i=l 



III. DETAILS OF THE SIMULATION 

We have performed Monte Carlo (MC) simulations 
at various /3-values between (3 = 2.2 and /3 = 2.55 
for various lattice sizes L. Consecutive configurations 
(considered to be statistically independent) were sep- 
arated by 100 sweeps, each sweep consisting of one 
local heatbath update followed by L/2 microcanoni- 
cal updates. In Table|T]we provide the full information 
about the field ensembles used in this investigation. 

For gauge fixing we employ the Z(2) flip opera- 
tion as discussed in For completeness we repeat 
the main information. The method consists in flip- 
ping all link variables U xll attached and orthogonal 
to a 3d plane by multiplying them with -1 £ Zi2). 
Such global flips are equivalent to non-periodic gauge 
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200 
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TABLE I: Values of (3, lattice sizes, number of measure- 
ments and number of gauge copies used throughout this 
paper. For the values of the lattice spacing see [34| (1 
GeV" 1 ~ 0.197 fm). 



transformations and represent an exact symmetry of 
the pure gauge action considered here. The Polyakov 
loops in the direction of the chosen links and averaged 
over the 3d plane obviously change their sign. There- 
fore, the flip operations combine the 2 4 distinct gauge 
orbits (or Polyakov loop sectors) of strictly periodic 
gauge transformations into one larger gauge orbit. 

The second ingredient is the simulated annealing 
(SA) method, which has been found computationally 
more efficient than the only use of standard overrelax- 
ation (OR) ESE1E1- The SA algorithm generates 
gauge transformations gix) by MC iterations with a 
statistical weight proportional to exp (41^ Fjj[g]/T) . 
The "temperature" T is an auxiliary parameter 
which is gradually decreased in order to maximize the 
gauge functional Fjj [g] . In the beginning, T has to 
be chosen sufficiently large in order to allow traversing 
the configuration space of gix) fields in large steps. 
As in Ref. El we have chosen T; n i t = 1.5. After each 
quasi-equilibrium sweep, including both heatbath and 
microcanonical updates, T has been decreased with 
equal step size. The final SA temperature has been 
fixed such that during the consecutively applied OR 
algorithm the violation of the transversality condition 



max 

x, a 



4 

E 



+A/2;m 



At 



(10) 



decreases in a more or less monotonous manner for 
the majority of gauge fixing trials until the condition 
(JTUJ) becomes satisfied with a unique e; or = 10~ 7 . A 
monotonous OR behavior is reasonably satisfied for a 
final lower SA temperature value Tg na i = 0.01 59]. 
The number of temperature steps has been chosen to 
be 1000 for the smaller lattice sizes and increased to 
2000 for the lattice size 30 4 and bigger. The final- 
izing OR algorithm using the standard Los Alamos 
type overrelaxation with the parameter value w = 1.7 
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requires typically a number of iterations varying from 
O(10 2 ) to O(10 3 ). 

In what follows we call the combined algorithm em- 
ploying SA (with finalizing OR) and Z(2) flips the 
'FSA' algorithm. By repeated starts of the FSA algo- 
rithm we search in each Z(2) Polyakov loop sector sev- 
eral times for the best ("6c ") copy. The total number 
of copies per configuration N copy generated for each 
/3-value and lattice size is indicated in Table HI In 
order to demonstrate the Gribov copy effect we can 
compare with the results obtained from the randomly 
chosen first ( "/c " ) copy. 

Some more details to speed up the gauge fixing pro- 
cedure are described in [421 ] . 



IV. INFLUENCE OF GRIBOV COPIES ON 
THE GLUON PROPAGATOR 

Our efficient gauge fixing procedure changes signif- 
icantly the momentum dependence of the gluon prop- 
agator D(p) in the IR- region in comparison with the 
result of the standard overrelaxation method. 

In Fig.[TJwe compare our be FSA results for the un- 
renormalized gluon propagator D(p) calculated on a 
44 4 lattice with those of the standard fc OR method 
obtained for an 80 4 lattice in Ref. [60]; all data pro- 
duced for (3 — 2.3 2 . One can see that our be FSA data 
points lie essentially below those of the fc OR method 
for momenta [p|^$0.7 GeV. Thus, we observe that the 
OR method with one gauge copy produces unreliable 
results for this range of momenta. Note that our lat- 
tice size is approximately twice as small as that in Ref. 

In the same figure we compare also with the fc SA 
results (no flips taken into account but much longer 
SA schedule applied) obtained for an 80 4 lattice in 
Ref. [6l[. The data look consistent for the momen- 
tum region |p|^0.3 GeV. At smaller momenta, e.g., 
at |p|^0.3 GeV the fc SA data points show the ten- 
dency to lie at somewhat lower values. The reason 
for this difference can be attributed to finite-volume 
effects (our lattice size is much smaller) or to uncer- 
tainties due to the lower statistics in [6l| . In any case 
we confirm that the Z(2) flips have the tendency to 
lower finite-size effects [571 ] - 

Let us define the Gribov copy sensitivity parameter 
A(p) as a normalized difference of the fc and be gluon 
propagators 



A(p) 



L» /C (p) - D bc (p) 
D bc (p) 



(11) 



2 Only momenta GeV are shown for data from |60l . note 

also that we are using a slightly different value for the lattice 
spacing than in |60( . 
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FIG. 1: Comparison of data obtained for be FSA gauge fix- 
ing with those obtained with the standard fc OR - method 
and the fc SA algorithm (all for p = 2.30). 
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FIG. 2: The momentum dependence of the Gribov copy 
sensitivity parameter A(p) for various lattices with physi- 
cal size aL ~ 3 fm. 



where the numerator has been obtained by averaging 
the differences between fc SA and be FSA propaga- 
tors calculated for every configuration and normalized 
with the be (averaged) propagator. 

In Figs. [31 El and 3] we show the momentum depen- 
dence of the Gribov copy sensivity parameter A(p) 
for different lattices with physical sizes aL ~ 3 fm, 
aL ~ 5 fm and aL ~ 7.3 fm, respectively. 

Evidently, the Gribov copy effect is rather strong 
in the deep IR-region. It does not disappear with 
rising [3, i.e. with the lattice spacing a becoming 
smaller. However, the effect decreases rapidly for ris- 
ing momentum. Moreover, it is encouraging to see 
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0.6 0.8 1.0 
Ipl; Gev 



method. Of course, we cannot exclude that other 
even more efficient gauge fixing methods can be in- 
vented which might bring us even closer to the global 
extremum of the gauge functional. But we are con- 
vinced that this will not change the conclusions given 
here. 

We would like to emphasize that the supposition to 
find the Gribov copies as close as possible to the global 
extremum of the Landau gauge functional differs sub- 
stantially from the recent claim to study those Gri- 
bov copies maximally enhancing the infrared asymp- 
totics of the ghost dressing function j45[. In as far 
these two different strategies really provide the differ- 
ent decoupling and scaling solutions, respectively, for 
the gluon and ghost propagators within the thermo- 
dynamic limit remains an interesting question. 



FIG. 3: The momentum dependence of A(p) for various 
lattices with physical size aL ~ 5 fm. 




Ipl; Gev 

FIG. 4: The momentum dependence of A(p) on a lattice 
with physical size aL ~ 7.3 fm. 



that the influence of Gribov copies for fixed momen- 
tum demonstrates clear tendency to decrease with in- 
creasing physical size aL. Thus, with our gauge fixing 
procedure we do not find any Gribov copy effect for 
momenta p^O. 6 GeV on lattices with aL ~ 3 fm, while 
on lattices with aL ~ 7.3 fm the gluon propagator for 
all momenta p^OA Gev is free of Gribov copy effect. 
This tendenc y is in accordance with a conjecture by 
Zwanziger i n 16211 and was seen already for smaller lat- 
tice sizes in [63j. 

Let us finally note that the above observations are 
valid for our FSA gauge fixing method (employed with 
rather large number of gauge copies), which is proved 
to be much more powerful than, e.g., standard OR 



V. NUMERICAL GLUON PROPAGATOR 
RESULTS 

In order to suppress lattice artifacts from the be- 
ginning we followed Ref. [3(| and selected the allowed 
lattice momenta as surviving the cylinder cut 

E^J(EM 2 <3. (12) 

Moreover, we have applied the "a-cut" (64[ fV — 
(2/a)a for every component, in order to keep close 
to a linear behavior of the lattice momenta p^ = 
(2nk li )/(aL), fc p G (-L/2.L/2]. We have chosen 
a = 0.5. Obviously, this cut influences large momenta 
only. 

We define the renormalized propagator D ren (p) 
according to the momentum subtraction schemes 
(MOM) by 

D ren (p, fi) = Z ren {n, l/a) D(p, l/a) (13) 
D ren {p = n) = 1/V 2 . (14) 

In practice we have fitted the bare propagators 
D{p, l/a) with an appropriate function (see Eq. (|15p 
below) and then used the fits for renormalizing D{p). 
But it has to be seen, that multiplicative renormal- 
izability really holds in the non-perturbative regime. 
For this it is sufficient to prove, that ratios of the 
renormalized (or unrenormalized) propagators ob- 
tained from different cutoff values l/a(/3) will not de- 
pend on p at least within a certain momentum in- 
terval Pmin,Pmax, where p m ax should be the maximal 
momentum surviving all the cuts applied. 

In what follows the subtraction momentum has al- 
ways been chosen to be fj, — 2.2 GeV. In Fig. [5] we 
show the momentum dependence of the renormal- 
ized gluon propagator D ren {p) at comparatively large 
momenta GeV) for five different lattice spac- 

ings but with (approximately) the same physical size 
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FIG. 5: The momentum dependence of the gluon propa- 



gator Dren(p) for five different lattice spacings and |p|;Cl 
GeV. The physical linear box size is aL ~ 3 fm. 
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FIG. 6: The same as in Fig.EJbut for |p|^l GeV. 



aL ~ 3 fm (for exact values see Table HJ. Evidently, 
in this momentum range the finite-spacing effects are 
rather small, at least for inverse bare coupling values 
A/gl = (3 > 2.2. 

In Fig. [6] for the same data set, we show the IR 
region only, whereas Fig. [7] presents the momentum 
dependence of the renormalized gluon dressing func- 
tion Z ren (p) — p 2 D ren (jp) for the same lattice spac- 
ings as in Figs. [5] and [6] In both cases we see that 
there are quite strong deviations at least for j3 = 2.2. 
The early turnover of the gluon propagator D ren (p) 
to an IR flattening is much less pronounced for higher 
(3 — values and thus, has to be attributed to lattice 
artifacts. We have to suspect that similar computa- 
tions done for SU(3) at the quite low value (3 — 5.7 
suffer from the same problem [43| • 

In Figs. [H [H and [10] we show an analogous set of 
data but now for the larger volume of about (5fm) . 
Also in this case we see that the continuum limit is 
fastly reached for [3 > 2.3 in the whole momentum 
range, whereas for (3 — 2.2 lattice artifact devia- 
tions occur, which become particularly strong in the 
infrared region |p|<l GeV. 

In order to estimate finite-volume effects we com- 
pare our be FSA data for the propagator D ren (j>) ob- 
tained for the same (3- values but different aL. In 
Fig. [11] we show the momentum dependence of the 
renormalized gluon propagators D ren (p) for two dif- 
ferent physical sizes aL at (3 = 2.4, and in Fig. [12] 
Dren (p) is presented for three different volumes at the 
somewhat stronger coupling (3 = 2.3. One can see that 
finite- volume effects are present only for the zero and 
minimal nonzero momenta and in the latter case they 
are rather small. Moreover, the IR flattening becomes 
visible only for the largest volume aL ~ 7 fm. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Ipl; Gev 

FIG. 7: The momentum dependence of the renormalized 
dressing function Z r en(p) for five different lattice spacings. 
The physical linear box size is aL ~ 3 fm. 



Evidently, our data for the zero-momentum propa- 
gator D ren (p = 0) show a clear tendency to decrease 
with increasing physical size aL. However, it remains 
difficult to speculate about the infinite volume limit 
of D ren (0). 

In the literature one can find quite a few functional 
forms suggested to describe the gluon propagator in 
the IR-region, most of them of purely phcnomenolog- 
ical origin, see e.g. [3(3]. For (3 > 2.3 we have fitted 
the momentum dependence of the gluon propagator 
with an ansatz describing a massive behavior in the 
infrared 

n ( \ — Ql i ° 2 i ° 3 



7 



1.5 




1.4 


h 


1.3 




1.2 
1.1 




< 1.0 


r ^ 


o 0.9 




0-8 




3 0.7 




|0.6 




0.5 




0.4 




0.3 




0.2 




0.1 




0.0 





I I I 111 

aL=5 fm 



A P=2.20; 24 4 
o |3=2.30; 30 4 
• |3=2.40; 42 4 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 
Ipl; Gev 



1.6 
1.4 
1.2 
S 1.0 



Q 

- a 0.8 



0.6 
0.4 
0.2 



ft 



A (3=2.20; 24 
o |3=2.30; 30 4 
• P=2.40; 42 4 



aL=5 fm 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 
Ipl; Gev 



FIG. 8: The momentum dependence of the renormalized 
gluon propagator D ren (p) for three different lattice spac- 
ings and |p|~l GeV. The physical linear box size is aL ~ 5 
fm. 



FIG. 10: The momentum dependence of the renormalized 
dressing function Z ren (p) for three different lattice spac- 
ings. The physical linear box size is aL ~ 5 fm. 
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FIG. 9: The same as in Fig.[8]but for |p|^l GeV. 
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FIG. 11: The momentum dependence of the renormalized 
gluon propagator D ren (p) at fixed /3 = 2.4 for two different 
physical volumes. 



cli, 0,2, a-3 and m are fit parameters. The fitting curves 
allow to compare quite easily the propagators ob- 
tained on different volumes. We show the fit results 
obtained for p>0.6 GeV in Figs. HUH and EES for 
the same three volumes as considered before. The fit 
parameters are provided in Table HT1 

In Fig. [TH] we compare the data together with the 
corresponding fit curves obtained at [3 = 2.3 for 
the three physical volumes. We see that in the range 
shown p > 0.6 GeV finite-size effects are small at 
least for physical linear lattice sizes aL > 5 fm. 
Finally we check the multiplicative renormalizability 



by presenting the data for the gluon propagator ratio 

R = D ren (p; (3; L)/D^(p; p = 2.55; L = 42) (16) 

in Fig. 1171 The relative deviations are below 3 percent 
such that we can say multiplicative renormalizability 
is safe for f3 > 2.4 and for the momentum range p > 
0.6 GeV. 
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FIG. 12: The momentum dependence of the renormalized FIG. 14: The momentum dependence of the renormalized 



gluon propagators D ren {p) at fixed /3 = 2.3 for three dif- 
ferent physical volumes. 



dressing function Z ren (p) as for Fig. [T3] but with linear 
box size aL ~ 5 fm. 




FIG. 13: The momentum dependence of the renormalized 
dressing function Z ren (p) fitted at the smallest available 
lattice spacing compared with the data of coarser lattices. 
The linear box size is aL ~ 3 fm. 



FIG. 15: The momentum dependence fit of the renormal- 
ized dressing function Z r an(p) for j3 — 2.3 and linear box 
size aL ~ 7 fm. 



VI. DEEP INFRARED REGION AND 
INFRARED MASS SCALE 

In the previous section we discussed our fits of the 
propagator in the range 0.6 ^p^S 3 GeV. An even 
more phcnomcnologically interesting range is the deep 
infrared region \p\ < 0.6 GeV. As it was discussed in 
Section [J the infrared mass scale is an important pa- 
rameter for phenomenological analyses. Lattice re- 
sults in principle can provide model independent in- 
formation on this subject. Our results obtained in this 
region for aL = 5 fm and j3 = 2.3 and 2.4 exhibit only 



a weak dependence on the lattice spacing (see Figs. 
[HI flU]) . Moreover, the finite-size dependence of D(0) 
seems to be moderate at f3 = 2.4 (compare with Figs. 
[TT1[T2|) . Therefore, let us speculate that our results for 
momenta below 0.6 GeV for j3 = 2.4 and aL — 5 fm 
are already close to the continuum and to the large- 
volume limit and that multiplicative renormalization 
can be assumed, too. Under these assumptions we 
have fitted our results at p^S 0.6 GeV separately. We 
used two fitting functions to fit our data in this range 
of momenta. The first one is just the pole-type prop- 
agator providing an effective gluon mass scale m g for 
p^0 
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FIG. 16: The renormalized dressing functions Z ren (jp) at 
[3 = 2.3. The data and the corresponding fits are shown 
for linear box sizes aL ~ 3, 5, 7 fm. 



FIG. 17: The momentum dependence of the i?-ratio for 
different j3 values and physical volumes. 



p 


L 


a\ 


a 2 


a3 


m 


X% 


2.20 


14 


0.62 


4.58 


49.5 


1.38 


5.0 


2.30 


18 


0.61 


4.03 


54.5 


1.37 


3.8 


2.40 


26 


0.56 


4.31 


50.8 


1.35 


1.5 


2.50 


36 


0.55 


4.30 


39.5 


1.30 


1.5 


2.55 


42 


0.56 


4.34 


45.0 


1.34 


1.2 


2.20 


24 


1.02 


1.33 


128.4 


1.53 


3.6 


2.30 


30 


0.59 


4.19 


61.0 


1.39 


2.6 


2.40 


42 


0.55 


4.39 


53.9 


1.37 


0.5 


2.30 


44 


0.67 


3.65 


72.2 


1.41 


1.6 



TABLE II: Values of the fit parameters in physical units 
(with dimension [a{\ = GeV°, [02] = GeV 2 , [03] = 
GeV 6 , [m] — GeV) and the corresponding Xdf ■ 



-Dpole(p) = 



A 



p 2 + m 2 



The other one is of the form of a Gaussian 



D 



gauss 



(17) 



(18) 



A, Trig, B,po, ttiib. are fit parameters. The results ob- 
tained for [3 — 2.4 and aL = 5 fm are compared in 

Fig. mi 

One can see that the pole-type momentum depen- 
dence given in Eq. (TP7)) is not suitable for fitting the 
gluon propagator in the considered range of momenta. 
This is supported by the large value of the xl n df 
parameter value which was found to be about 14 if 




MC; |3=2.40; 42 
Pole fit 
Gauss fit 



5.4 

5.2 

5.0 

4.8 

„ 4.6 

u 4.4 
O 

.. 4.2 

@ 4.0 
Q S 3.8 
3.6 
3.4 
3.2 

3.0^ 

0.2 0.3 0.4 

Ipl; Gev 

FIG. 18: The renormalized gluon propagator at f3 = 2.4 
on aL = 5 fm lattice together with the infrared fits acc. 
to Eqs. (|17I18|) . 



p = is included and about 7 if this point is ex- 
cluded. Therefore, the definition of an effective gluon 
mass m g via the tree-level Eq. (jTTJ) turns out to be 
problematic. 

On the contrary, the Gaussian-type momentum de- 
pendence given in Eq. (fT8|) nicely fits the data for 
the gluon propagator in the range p ^$ 0.6 GeV with 
X 2 /ndf about 0.7 (0.5) with p — included (ex- 
cluded). Therefore, the Gaussian functional form al- 
lows to define an infrared (mass) scale roiR. The nu- 
merical fit result is toir = 0.69(3) MeV forp = taken 
into account in the fit and 0.68(4) MeV, if p = is 
ignored. 



10 



It is also interesting to note that our data suggest 
the existence of a maximum of the gluon propagator at 
a non-zero momentum p . We have already observed 
this for SU (2) in four dimensions in our previous work 
[42| , where we studied the propagator on large lattices 
at (3 = 2.2, i.e. rather far from the scaling region. 
Our new results reported for (3 — 2.3 and 2.4 (see 
Figs. I11|T2"|) show that the maximum seems to persist 
in the continuum limit. With our fit function Eq. (fTS)) 
we obtained its position at pa — 180(15) MeV. The 
mere existence of the maximum already contradicts a 
simple effective gluon mass prescription ~ 1/ (j> 2 + m 2 g . 



VII. CONCLUSIONS 

In this work we investigated numerically the renor- 
malized Landau gauge gluon propagator D ren {p) in 
the pure gauge SU(2) lattice theory. The main goal 
of this study was to study the approach to the con- 
tinuum limit, especially in the infrared region |p|P$l 
GeV. 

In order to disentangle finite-spacing from finite- 
volume effects we calculated the propagators on lat- 
tices with physical size aL equal approximately 3 fm 
at various /3-values in the range of j3 = 2.2, . . . , 2.55 
and on lattices with aL w 5 fm for (3 = 2.2, . . . , 2.4 
(see Table Q] for details) . Calculations were made also 
on aL w 7.3fm lattices at (3 = 2.3. Our lattice vol- 
umes varied from L 4 = 14 4 to 44 4 . For physical vol- 
umes (aL) 4 ~ (3fm) 4 and (5fm) 4 , we have checked 
the scaling behaviour assuming the lattice spacing a 
to depend on (3 as determined from the string tension. 
The comparison of the renormalized propagators cal- 
culated for different physical volumes then allowed to 
estimate the influence of the finite (physical) volume 
in the infrared regime. 

Special attention has been paid to the dependence 
on the choice of Gribov copies. In our previous papers 
[57L Ir35| we have seen that the finite- volume behavior 
of the gluon propagator (and not only of the ghost 
propagator) is sensitive to the way how the Landau 
gauge is fixed. We found indications that by enlarg- 
ing the gauge orbits by Z(2)-flip operations and by 
applying the simulated annealing method with con- 
secutive overrelaxation ('FSA' algorithm) the volume 
dependence becomes suppressed. However, our former 
investigations went immediately to largest accessible 
lattice volumes by employing coarse lattices {(3 = 2.2), 
such that the continuum limit remained to be studied. 

Our findings can be summarized as follows. 

1) In the region |p|^0.6 GeV we observe very nice 
agreement between renormalized propagators D ren (p) 
obtained for different lattices with (3 > 2.4. For larger 
momenta scaling holds for even smaller /3 values. In 
contrast, in the deep infrared region (|p|<0.6 GeV) 
the scaling violation is quite strong, especially for [3 — 



2.2. However, with increasing j3 finite-spacing effects 
rapidly decrease: our data seem to 'converge' to some 
limiting curve of the renormalized propagator D ren (p) 
(compare, e.g., propagators at (3 = 2.2 and j3 = 2.55), 
thus indicating the approach to the continuum limit 
in the fixed physical volume (aL) 4 . 

2) Using our gauge fixing procedure we observed finite 
(physical) volume effects only for zero and minimal 
nonzero momenta. We made this observation for (3 = 
2.2 in [42j], here we confirm it for (3 — 2.3 and 2.4 and 
thus it can be extended to the continuum limit. 

3) In our previous papers [!2,[52l we calculated gluon 
propagators on various lattices at (3 = 2.2. We ob- 
served the appearance of a maximum at a non-zero 
value of the momentum p on lattices with compar- 
atively large volume (aL > 6.7 fm). Our new data 
obtained at larger /3-values do confirm this observa- 
tion. Fitting the propagator calculated at j3 = 2.4 
by the fitting function Eq. I|18p we obtained our best 
estimation for the position of this maximum as po = 
180(15) MeV. This number might slightly change in 
the continuum limit. Also we observed that the zero- 
momentum gluon propagator D(0) has a tendency to 
decrease with growing lattice size L. We did not try 
here to extrapolate its infinite- volume value. 

4) The effective gluon mass m g has been employed as 
an important parameter in various phenomenological 
analyses. We fitted our results at p^0.6 GeV using 
our data for the gluon propagator at (3 — 2.4 on lat- 
tices with aL ps 5. It was found that the pole- type 
momentum dependence given by Eq. (fl7)) does not 
provide an adequate description of our data at small 
momenta, which makes it problematic to define this 
parameter from the Landau gauge gluon propagator. 
On the contrary, a Gaussian-type behavior given by 
Eq. (fTB|) fits the data nicely. It allows to define an 
alternative infrared (mass) scale miR = 0.69(3) GeV 
describing the approach to the infrared limit. Its con- 
sequences for phenomenological applications remain 
to be seen. 

5) We confirm that the Gribov copy influence is very 
strong in the deep infrared region. Comparing our 
be FSA results calculated on a 44 4 lattice with those 
of the standard fc OR method obtained for an 80 4 
lattice in Ref. [6(| (see Fig. [T|) we found that the 
OR method with one gauge copy produces unreliable 
results for momenta |p|^;0.7 GeV. We conclude that 
fc OR method should be applied exclusively to large- 
momentum studies. 

Our FSA method provides systematically higher 
values of the gauge fixing functional as compared to 
the standard OR procedure. We studied the Gribov 
copy effect for this method as well, generating up to 
80 gauge copies for every configuration. We found for 
fixed physical volume the Gribov copy sensivity pa- 
rameter A(p) only weakly to depend on the lattice 
spacing a. Therefore, the quality of the gauge fixing 
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procedure in the study of gauge dependent observables 
remains important, at least, in the deep infrared. 

At the same time the influence of Gribov copies 
demonstrates clear tendency to decrease for fixed mo- 
mentum with increasing physical size aL. This ten- 
dency is in accordance with a conjecture by Zwanziger 
in 163 and was seen already for smaller lattice sizes 
in 16311. 
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